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A NOISE SUPPRESSOR 



FIELD OF THE INVENTION 

5 This invention relates to noise suppression and is particularly, but not exclusively, 
related to noise suppression in a speech signal picked up by a mobile terminal 
such as a mobile phone. 

BACKGROUND OF THE INVENTION 

10 

When a communications terminal is used to make a record of or to transmit a 
speech signal containing speech, it is inevitable that its microphone will pick up 
environmental or background noise from the environment in which a speaking 
person is located. The background noise reduces the ability of a listener to hear or 

15 understand the speech and in some cases, if the noise level is sufficiently high, 
prevents the listener from hearing anything other than the background noise. In 
addition, such background noise may have a negative effect on the performance 
of digital signal processing systems in the communications terminal or in an 
associated communications network, such as speech coding or speech 

20 recognition. Typically, noise suppression systems are incorporated in 
communications terminals and communications networks to limit the effect of 
background noise. 

Noise suppression has been well known for a number of years. Many different 
25 approaches and methods have been proposed to achieve three main ends: 

(i) suppressing the noise significantly while preserving good speech quality; 

(ii) rapid convergence to the optimal solution independent of the nature of the 
processed noise; and 

(iii) improving speech intelligibility for very low speech-to-noise (SNR) ratios. 
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One noise suppression method based on the linear Minimum Mean Squared Error 
(MMSE) criteria will be described with reference to Figure 1 . The method operates 
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on a noisy speech signal x(t) containing a speech signal s(t) and a noise signal 
n(t) such that x(t) = s(t) + n(t) . The noisy speech signal x(t) is in the time domain. 
It is converted into a sequence of frames having consecutive frame numbers k 
using a windowing function. The frames are then each transformed into the 
5 frequency domain using a Fast Fourier Transform (FFT) in block 10 so as to 
produce a sequence of noisy speech frames where noisy speech signal X(f,k) in 
the frequency domain contains a speech signal S(f,k) and a noise signal N(f,k) 
such that X(f,k) = S(f,k) + N(f,k). The frames in the frequency domain 
comprise a number of frequency bins / . In the frequency domain, the MMSE 
10 approach involves minimising the following error function: 

e\f,k) = E{s(f,k)-S(f,k)y(s(f 9 k)-S(f,k))*} (1) 

where E{} is the expectation operator, (*) denotes complex conjugation and 
15 S(f,k) represents a linear estimate of the input speech signal. The error 
e 2 (f,k) defined by Equation 1 represents the squared difference between the true 
speech component contained within the noisy speech signal and the estimate of 
that speech component, S(f,k), i.e. the estimate of the noise-free speech 
component. Thus, minimisation of £ 2 (f,k)\s equivalent to obtaining the best 
20 possible estimate of the speech component. S(f,k) is given by: 

S(/,*) = G(/,*) ■*(/,*) (2) 

where G(f,k) is a gain coefficient. The corresponding solution of the minimisation 
25 of e 2 (f,k) for each frame takes the form of a computation of the gain coefficient 
G(f,k) which is multiplied by the associated input frequency bin of that frame to 

produce the estimated noise-free speech component S(f,k). This gain 

coefficient, known as the frequency domain Wiener filter, is given by the ratio 
below: 
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G(f,k) = 



E{S(f,k)X\f,k)} 
E{x(f,k)X*(f,k)} 



(3) 



The Wiener filter G{f,k) , is generated for each frequency bin / of each frame. 

5 

The noise-suppressed frames are then transformed back into the time domain in 
block 14 and then combined together to provide a noise suppressed speech signal 
5(0 . Ideally, s(t) = s(t) . 

10 When deriving the Wiener filter, the MMSE approach is equivalent to the 
orthogonality principle. This principle stipulates that, for each frequency, the input 
signal X(f,k) is orthogonal to the error S(f,k)-S(f,k) . This means that: 



Because the estimation process is linear, by estimating the signal component of a 
noisy signal that contains a signal component and a noise component, an estimate 
of the noise N(f,k) is also effectively obtained. Furthermore, the following 
orthogonality relationship will also be true: 



where N(f,k) indicates the noise estimate. It also follows that for every 
frequency, the following equality applies: 



e{(s(/,*)-S(/,*))- **(/,*)}= 0 



(4) 
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E{(N(f,k)-N(f,k))x\f,k)}=0 



(5) 
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S(f,k) - S(f,k) = N(f,k) - N(f,k) 



(6) 
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that is, the error associated with the estimate of the noise component N(f,k) is the 
same as the error associated with the estimated noise-free speech component 



In the remainder of this document, the following notation will be adopted: P uv (f,k) 
is the cross power spectral density between U(f,k) and V(f,k) 
(Pvv(f,k) = Ep(f 9 k)-V*(f 9 k)}).P uu (f,k) is the power spectral density (psd) of 

U(f,k) (P uu (f,k) = E{U(f,k)-U\f,kj\). 

As a consequence of the above-mentioned orthogonality principle, it is possible to 
derive an expression for the cross psd P sx (f,k), required in order to compute the 
Wiener filter described by Equation 3: 

P sx (f,k) = E{{x(f,k)-N(f,k)).X\f,k)} (7) 
Moreover, the cross psd P^ (f,k) is given by: 

P m (f,k) = E{(X(/,*) - S(/,*)). X * (/,*)} (8) 

Having in mind the trivial equality P xx (f,k) = P sx (f,k) + P„ x (f,k) i Equations 3, 6, 
7 and 8 introduce and illustrate an idea of adaptive calculation since the Wiener 
filter (P sx (f \k)) in Equation 3 depends on the estimated signal S(f,k) 
(6,7) and (8). 

When a minimum is reached, the expression describing the error in Equation 2 
takes the following form: 



S(f,k). 



*L </.*) = 



O) 
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It is evident that minimum error, that is (/,*), is e q ual to z oro onl Y if the 
desired signal S(f,k) is completely coherent with the input signal X(f,k) (that is, 
P m (f,k) tends to zero). This is desirable. Otherwise, there is an error when 
5 applying the Wiener filter. The upper limit of this error is P ss (f 9 k). This is 

undesirable. In other words, an error-free result can only be obtained if there is 
actually no noise in the input signal X(f,k) . For any finite noise level, a finite error 
is obtained. It follows that the worst case error occurs when there is no speech 
signal S(/,*)in X(f 9 k). 
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* SUMMARY OF THE INVENTION 



According to a first aspect of the invention there is provided a method of 
suppressing noise in a signal containing noise to provide a noise suppressed 
15 signal in which an estimate is made of the noise and an estimate is made of 
speech together with some noise. 



Preferably the signal comprises speech. 



20 Preferably the level of the noise included in the estimate of the speech together 
with some noise is variable so as to include a desired amount of noise in the 
noise-suppressed signal. 

Preferably the level of the noise provides an acceptable level of context 
25 information. 



Preferably the level of the noise is below the mask limit of the speech and so is not 
audible to a listener. Alternatively the level of noise approaches the mask limit of 
the speech and so some noise context information is left in the signal. 
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Preferably the method does not suppress noise if the signal to noise ratio is 
sufficiently high so that the level of noise already provides an acceptable level of 
context information or is already below the mask limit. 

5 Preferably the estimated noise is power spectral density. 

According to a second aspect of the invention there is provided a method of 
producing a gain coefficient for noise suppression in which a first estimation of the 
gain coefficient is made adaptively and this first estimation is used to produce a 
10 noise estimation which is then used to produce a second estimation of the gain 
function. 

In this respect, the invention provides an important advantage. It effectively 
eliminates the need for a Voice Activity Detector (VAD) in a noise suppressor 

15 implemented according to the invention. A VAD is basically an energy detector. It 
receives a noisy speech signal, compares the energy of the filtered signal with a 
predetermined threshold and indicates that speech is present in the received 
signal whenever the threshold is exceeded. In many speech encoding/decoding 
systems, particularly in the field of mobile telecommunications, operation of the 

20 VAD changes the way in which background noise in a speech signal is processed. 
Specifically, during periods when no speech is detected, transmission may be cut 
and so-called "comfort noise" generated at the receiving terminal. Thus use of 
such discontinuous transmission and voice activity detection schemes may 
complicate the use of noise suppression and lead to unwanted effects. Elimination 

25 of the need for a voice activity detector and the creation of a noise suppression 
scheme that automatically adapts to changes in noise conditions is therefore 
highly desirable. Because the invention introduces a method of noise suppression 
in which an estimate of both speech and background noise is obtained, there is 
effectively no need to make a decision as to whether an input signal contains 

30 speech and noise or just noise. As a result the VAD function becomes redundant. 

Preferably the first estimation is used to up-date the estimated noise. 
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According to other aspects of the invention, there is provided a noise suppressor 
operating according to the first aspect of the invention, a noise suppressor 
operating according to the second aspect of the invention, a noise suppressor 
5 operating according to the first and the second aspects of the invention, a 
communications terminal comprising a noise suppressor according to the first 
and/or second aspects of the invention and a communications network comprising 
a noise suppressor according to the first and/or second aspects of the invention. 

10 Preferably the communications terminal is mobile. Alternatively, the invention may 
be used in a network or fixed communications terminal. 



According to another aspect of the invention there is provided a method of 
calculating a Wiener filter in which an estimate is made of speech and background 
15 noise and the noise is far enough below the speech so that it is wholly or partially 
masked below the audible level or perception of a user. 

Preferably the method is for noise suppression in the frequency domain. It may 
comprise calculating the numerator and denominator of a Wiener filter to be used 
20 for a noise reduction system. The noise suppression system described in this 
document is particularly suitable for application in a system comprising a single 
sensor such as a microphone. 

Preferably the filter is a Wiener Filter. Preferably it is based on an estimate of a 
25 periodogram comprising a combination of speech and noise. Preferably the 
method involves continuous up-dating of noise psd. 



BRIEF DESCRIPTION OF THE DRAWINGS 
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An embodiment of the invention will now be described by way of example only with 

reference to the accompanying drawings in which: 

Figure 1 shows a mobile terminal according to the invention; 
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Figure 2 shows a noise suppressor according to the invention; 

Figure 3 shows the frequency and sound level dependent masking effect of the 

human auditory system 

Figure 4 shows a block diagram of an algorithm according to the invention; and 
5 Figure 5 shows a functional block diagram of an algorithm according to the 
invention. 

DETAILED DESCRIPTION 

10 In the following the symbol P generally represents power. Where it is primed, that 
is P\ it represents a periodogram and where it is not primed, that is P, it 
represents a power spectral density (psd). In accordance with their generally 
accepted meanings, the term "periodogram" is used to denote an average 
calculated over a short period and the term power spectral density is used to 

1 5 represent a longer term average. 

An embodiment of a mobile terminal 10 comprising a noise suppressor 20 
according to the invention will now be described with reference to Figure 1. Figure 
1 corresponds to an arrangement of a mobile terminal according to the prior art 
20 although such prior art terminals comprise conventional prior art noise 
suppressors. The mobile terminal and the wireless communications system with 
which it communicates operate according to the Global System for Mobile 
telecommunications (GSM) standard. 

25 The mobile terminal 10 comprises a transmitting (speech encoding) branch 12 and 
a receiving (speech decoding) branch 14. In the transmitting (speech encoding) 
branch 12, a speech signal is picked up by a microphone 16 and sampled by an 
analogue-to-digital (A/D) converter 18 and noise suppressed in the noise 
suppressor 20 to produce an enhanced signal. This requires the spectrum of the 

30 background noise to be estimated so that background noise in the sampled signal 
can be suppressed. A typical noise suppressor operates in the frequency domain. 
The time domain signal is first transformed into the frequency domain which can 
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be carried out efficiently using a Fast Fourier Transform (FFT). In the frequency 
domain, voice activity is distinguished from background noise and when there is 
no voice activity, the spectrum of the background noise is estimated. Noise 
suppression gain coefficients are then calculated on the basis of the current input 
5 signal spectrum and the background noise estimate. Finally, the signal is 
transformed back to the time domain using an inverse FFT (IFFT). 

The enhanced (noise suppressed) signal is encoded by a speech encoder 22 to 
extract a set of speech parameters which are then channel encoded in a channel 
10 encoder 24, where redundancy is added to the encoded speech signal in order to 
provide some degree of error protection. The resultant signal is then up-converted 
into a radio frequency (RF) signal and transmitted by a transmitting/receiving unit 
26. The transmitting/receiving unit 26 comprises a duplex filter (not shown) 
connected to an antenna to enable both transmission and reception to occur. 

15 

A noise suppressor suitable for use in the mobile terminal of Figure 1 is described 
in published document W097/221 1 6. 

In order to lengthen battery life, different kinds of input signal-dependent low 
20 power operation modes are typically applied in mobile telecommunication 
systems. These arrangements are commonly referred to as discontinuous 
transmission (DTX). The basic idea in DTX is to discontinue the speech 
encoding/decoding process in non-speech periods. Typically, some kind of comfort 
noise signal, intended to resemble the background noise at the transmitting end, is 
25 produced as a replacement for actual background noise. 

The speech encoder 22 is connected to a transmission (TX) DTX handler 28. The 
TX DTX handler 28 receives an input from a voice activity detector (VAD) 30 which 
indicates whether there is a voice component in the noise suppressed signal 
30 provided as the output of noise suppressor block 20. If speech is detected in a 
signal, its transmission continues. If speech is not detected, transmission of the 
noise suppressed signal is stopped until speech is detected again. 
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In the receiving (speech decoding) branch 14 of the mobile terminal, an RF signal 
is received by the transmitting/receiving unit 26 and down-converted from RF to 
base-band signal. The base-band signal is channel decoded by a channel decoder 
5 32. If the channel decoder detects speech in the channel decoded signal, the 
signal is speech decoded by a speech decoder 34. 

The mobile terminal also comprises a bad frame handling unit 38 to handle bad, 
that is corrupted, frames. 



The signal produced by the speech decoder, whether decoded speech, comfort 
noise or repeated and attenuated frames is converted from digital to analogue 
form by a digital-to-analogue converter 40 and then played through a speaker or 
earpiece 42, for example to a listener. 



Further details of the noise suppressor 20 are shown in Figure 2. It comprises a 
Fast Fourier Transform, a gain coefficient or Wiener filter calculation block and an 
Inverse Fast Fourier Transform. Noise suppression is carried out in the frequency 
domain by multiplying frames by gain coefficients/Wiener filters. 



The operation of the noise suppressor 20 will now be described. According to the 
invention, rather than attempting to estimate the "true" speech component 
S(f,k) in a noisy speech signal, a Wiener filter is used to estimate a combination 
of speech and a certain amount of noise according to the relationship 
25 S(f,k) + ^-N(f,k) . The modified Wiener filter thus created takes the form: 
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G(f,k) = 




(10) 



P sx {f,k) + Z-P m (f,k) 
P sx {f,k) + P m {f,k) 
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Assuming that the speech and noise component are uncorrected (that is, the 
cross psd between the speech and noise components must be equal to zero, 
p sN(f> k ) = °)> Equation 10 can be re-expressed in the form: 



5 G{ f 9k)= P Ss(f>V + Z' P NN(f>k) 

P ss {f,k) + P NN {f,k) 



The role of the factor £ is explained below. 



As explained earlier, the main advantage of estimating a combination of speech 
% 10 and a certain amount of noise is that there should be less error associated with the 
H estimation. This benefit becomes further apparent in connection with Equation 12, 

SBfSSl 

[y presented below, which defines the minimum error obtained in this situation: 
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It can now be understood that as P NN (f,k) tends to zero, equation 12 tends to 

zero and so the error tends to zero as in the case of the prior art. In common with 
the prior art, this is desirable. However, since Equation 12 includes the factor of 

(l-£) 2 it reaches zero more quickly than in the case of the prior art. On the other 
20 hand, as P NN (f,k) increases, e^ n tends to (l-£) 2 • P ss (f,k). In common with the 

prior art, this is undesirable. However, the error provided by the method according 
to the invention is always smaller than that provided by the prior art method 
described earlier. This advantage arises because the multiplying factor (l-£) 2 
always serves to reduce the amount of error. Furthermore, the factor (l-£) 2 can 
25 be minimised by setting £ to an appropriate value, in which case the error is 
further minimised. 
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In the invention it has been recognised that the value of £ can be determined to 
achieve the following results: 

1. To provide a value of the product £-P NN (f,k) which is "masked" by 
5 p ss(f> k )- Ev© n though an estimate of combined speech and noise is 

computed, a listener will hear only speech because the product £-P NN (f,k) 

will be below his audible level of perception. In this way, advantage is taken 
of the properties of the human auditory system, allowing the speech 
periodogram to be calculated together with the maximum of masked noise 
10 periodogram. When £ is being applied to achieve this result, it is referred to 

as . 

The "masking" effect is a property of the human auditory system which 
effectively sets a frequency dependent and sound level dependent lower 

15 limit or threshold on auditory perception. Thus, any noise or speech 

components below the masking threshold will not be perceived (heard) by 
the listener. It is generally accepted that the masking threshold is 
approximately 13dB below the current input level, irrespective of frequency. 
This is illustrated in Figure 3. According to the invention, in order to 

20 estimate the pure speech signal (that is, when trying to eliminate all the 

background noise), it is sufficient to estimate the pure speech signal 
together with that part of the noise just below the masking threshold. 

2. To allow the level for noise reduction at the output to be freely chosen. This 
25 can be used to restore near-end context to the signal for the far-end 

listener. When £ is being applied to achieve this result, it is referred to as 

£ 2 . This means that £ may be chosen in such a way as to ensure adequate 
noise suppression, but also to permit a certain noise component to remain 
in the signal at the receiving terminal, such that the background noise 
30 appears to naturally represent the background noise present in the 

environment of a transmitting terminal. In other words it is possible to 
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choose a value of £ such that the noise component in a noisy speech 
signal is not completely eliminated due to the masking effect. 



In practical situations, speech signals are non-stationary and therefore require 
5 short-term estimation. Thus, instead of using psd functions, as shown in Equation 
11, certain terms are replaced with periodograms. Noise may be also non- 
stationary, but it is generally considered to be stationary, so long-term estimation 
may be still be used. Hence, the form of the desired Wiener filter is: 



10 G(f,k)= P s s y ,k) + ^ P ™ (f ' k) (13) 
P' ss (f,k) + P NN {f,k) 



It should be noted that it is also possible to use the background noise power 
spectral density term P NN (f,k) in the denominator of Equation 13. It should also be 

appreciated that when £ = is used in Equation 13 above, the term 
15 P f ss {f ,k) + ^ x *P' NN (f,k) represents a combination of the speech periodogram and 
the masked noise periodogram and when £ = ' s used, the term 
p ss(f + ' p nn(/^^ represents a combination of the speech periodogram and 
the permitted noise periodogram. The denominator P' ss {f ,k) + P NN (f ,k) is 
composed of the speech periodogram and the noise psd, respectively. 

20 

Calculation of the Wiener filter for a current frame k is based on a previous frame 
k-\ as follows. The noise psd P NN (f,k-l), the speech periodogram P' ss {f,k-\) 

and the number of frames T(f,k-1) for time averaging of previous frames are 
known. For the current frame k , a combination of the input speech and the noise 
25 periodogram \X(f,k)\ 2 is also known. Rather than P NN (f 9 k-l), R NN (f,k-\) or 
L NN (f,k-l) may be used if square root or logarithmic measures are employed, as 
described later in this description. 
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An eight-step algorithm is used to calculate the Wiener filter. The eight steps are 
shown in Figure 4 and are described below. 

Stepl: Estimation of a combination of the speech and the noise 
5 periodogram P s ' s (f,k) 

This periodogram is calculated as follows: 

p; s (/, k) = a • P; s (/, k - 1) + (1 - a) • |X (/, kf (14) 

10 

It should be noted that P^ s {f,k) is based on the previous periodogram of speech 

Pss(f> k -1) anc * an amount of the current noisy speech signal \X(f,k)\ 2 , 
determined by a factor a . The value of a is chosen to provide the greatest 
possible contribution from the current speech component |s(/,A;)| 2 of the noisy 

15 speech SIGNAL \X(f,k)\ 2 , but it is limited to ensure that the factor 
(l-a)-\N(f,k)\ 2 , which represents the amount of the current noise signal that will 

be included, is masked by the sum a- Ps S (f,k-l) + (l-a)-\S(f,k)\ 2 which 

represents an estimate of the current speech periodogram. Therefore, it should be 
appreciated that it is necessary to re-calculate the forgetting factor a for every 
20 frequency bin / of every frame k. It should also be noted that the factor (I -a) 
referred to in Equation 14 is analogous to 

Practically, step 1 is implemented by first estimating the current speech 
periodogram using the spectral subtraction method described in "Suppression of 
25 Acoustic Noise in Speech Using Spectral Subtraction", IEEE Trans. On Acoustics 
Speech and Signal Processing, vol. 27, no. 2, pp. 113-120, April 1979. Then the 
masking level is set at a value which is approximately 13dB below the estimated 
speech periodogram level. The noise periodogram is estimated in same way as 
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the speech periodogram. The value of a is then computed using the mask, the 
noise periodogram and the input periodogram. 

Step 2: Estimation of a combination of speech and noise psd P^ (f,k) 
This psd represents the total power of the input and is estimated by: 



P xx (f,k) = a* 



P; s (f,k-l) + -P NN (f,k-l) 



a 



(l-a)-\X(f,kf 



(15) 



10 This psd combines short term averaging (a periodogram for speech) together with 
long term averaging (a psd for noise). 



Step 3: Estimation of the Wiener Filter 
s 15 The Wiener filter of Equation 1 1 can be re-written in the following form: 



G l (f,k)= I J s{f ' k) (16) 



and so can be calculated from the results of Equations 14 and 15. Since 
20 S l (f,k) = G ] (f,k)-X(f,k) J it should be understood that the estimated speech 
contains the speech and the masked part of the noise. The minimum value 
for the gain G x {f,k) is set to (1-or). 

Step 4: Updating of the noise psd P NN (f,k) 

25 

To update the noise psd, the theoretical result presented in Equation 8 is used, 
replacing the product (X(f,k)-S(f,k))-X*(f,k) with the product 

(l-G x (f 9 ky)-\X(f,k)\ 2 where necessary. The following three methods can be 
used: 




an 
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(i) power psd estimation; 

(ii) square root psd estimation; and 

(iii) logarithm psd estimation. 

5 In all of the methods described below, X represents a forgetting factor between 0 
and 1 . 

(i) Power psd estimation 

10 This method uses the orthogonality principle and is based on the Welch method 
described in "The Use of Fast Fourier Transform for the Estimation of Power 
Spectra: A Method Based on Time Averaging Over Short, Modified 
Periodograms", IEEE Trans. On Audio and Electroacoustics, vol. AU-15, n. 2, pp. 
70-73, June 1967. It uses a technique known as "exponential time averaging", 

1 5 according to which: 

P NN (f,k) = X*P NN (f,k-l) + (l-X)*(l-G l (f,k))-\X(f,kf (17) 
where G x (f,k) is the Wiener filter calculated according to equation 16. 

20 

(ii) Square Root psd estimation 

This method uses a modification of the Welch method and is based on amplitude 
averaging: 



25 



R NN (/, k) = X ■ R m (f 9 k - 1) + (1 - X) • V(l - G x (/,*)) • \X (/, *)| 
(f,k) = R NN (f,k)-R NN (f,k) 



R m (f,k) represents an average noise amplitude. 
30 (iii) Logarithmic psd estimation 
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This method uses time averaging in the logarithm domain: 

^L NN (f,k) = A-L NN (f,k-l) + (l-Z)Log[(l-G 1 {f,k))\X(f,k)\ 2 ] 
(/,*) = exp[L„„ (/,*) + /] 

L NN (f,k) refers to an average in the logarithmic power domain, y is Euler's 
constant and has a value of 0.5772156649. 



In each of the three methods described above, the forgetting factor X plays an 
10 important role in the updating of the noise psd and is defined to provide a good 
psd estimation when noise amplitude is varying rapidly. This is done by relating X 

to differences between the current input periodogram \X(f,k)\ 2 and the noise psd 

P NN (f,k-l) in the previous frame. X depends on a value T(f,k) which defines 

the number of frames used for time averaging and is determined as follows: 

15 

if \X(f 9 kf >\0P m (f,k- 1) T(f 9 k) = 5 
<elseif\X(f,k)\ 2 <0.1-P NN (f,k-l) T(f 9 k) = 5 (20) 
else T(f 9 k) = Min[T(f,k - 1) + 1, 20] 



and X is derived from T(f 9 k) as follows: 

X= T(f,k) (21) 

20 

It should be noted that it is necessary to re-calculate the forgetting factor X for 
each frame k and for every frequency bin / . Clearly, as X is required in step 2, it 
needs to be calculated so that it is available for that step. It should also be 
appreciated that because the noise psd is updated continuously, this removes the 
25 need to have a voice activity detector in the noise suppressor 20. 



Step 5: 



Estimation of Current Speech Periodogram P' ss {f,k) 



+ 



The current speech periodogram P' ss {f,k) plays an important role in the algorithm. 
It is estimated for a current frame so that it can be used in a next frame, that is in 
Equations 14 and 15. As explained below, P/ 5 (f,k) should only contain speech 
5 and should not contain any noise. 

Effectively, after obtaining an estimate of speech amplitude S(f,k) in step 3, this 
step requires estimation of Ps S (f,k) which represents the current speech 
periodogram. 



It is widely accepted that P' ss {f,k) can simply be replaced with the squared 



Unfortunately, a good estimate S(f,k) does not actually imply that a good 
estimate for \S(f,k)\ 2 can be obtained by simply taking the square. Thus, the 
15 method according to the invention seeks to obtain a more accurate estimate 
Pss(f> k ) of \S(f>k)\ 2 by applying the MMSE criterion. 

Examining the combined speech and noise periodogram, it can be seen that: 

20 Y(f,k) = \X(f,kf = \S(f,k)\ 2 +\N(f,kf+S*(f,k)-N(f,k) + S(f,k).N*(f,k). 

Thus a good estimate of \S(f,k)\ 2 may be obtained by minimising the following 
error (MMSE criterion): 
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estimated speech amplitude, that is: P' ss (f,k) = \s(f 9 kj{ estimate of \S(f,k)\ 2 . 



25 




(22) 



where H(f,k) *\X(f,k)\ 2 represents an estimate of the speech periodogram 

\S(f,kf. 
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Direct solution of Equation 22 requires solution of higher order equations, but the 
solution can be simplified by assuming that the speech and noise are Gaussian 
processes, uncorrelated with zero means, to provide an approximation of the 
corresponding Higher Order Wiener filter H(f,k) . The approximation used in this 

method is presented in Equation 23 below. (It should be appreciated that different 
approximations may be used at this stage without departing from the essential 
features of the inventive principle). 



10 H(f k)= 3 SNR(f,k)*SNR(f,k) + SNR(f,k) 

3-SNR(f,k)-SNR(f,k) + 6-SNR(f,k) + 3 



Here, SNR(f,k) refers to the signal-to-noise ratio and is calculated as follows: 

SNR(f,k)= G ]?'*] (24) 

1-Ctj (J ,k) 



15 Equation 24 is the reciprocal of a well-known function relating the Wiener filter and 
the signal-to-noise ratio. (Wiener = SNR/(SNR+1)) 

Consequently, the speech periodogram is calculated as follows: 
20 /£(/,*) = H(f,k)-\X{f,kf (25) 
Step 6: The Amplification Function 

In conditions of high SNR, when the speech component of the noisy input signal is 
25 large compared with the noise component, the estimated Wiener filter G } (f,k) 

tends to 1. Furthermore, when the speech to noise ratio is high, G x {f,k) can be 
estimated comparatively accurately. Thus, there is a good degree of certainty that 
the Wiener filter determined in Step 3, offers optimal filtering and provides an 
output containing a highly accurate estimate of the speech S x (f) with a residual 



# 
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amount of (masked) noise. As the gain of the filter is close to 1 in this situation, it is 
advantageous to provide a small amount amplification to bring the gain still closer 
to 1. However, the additional amplification should also be limited to ensure that 
Wiener filter gain does not exceed 1 in any circumstance. 

5 

On the other hand in conditions where the speech component in the noisy input 
signal is small compared with the noise component, the opposite is true. The 
Wiener filter gain is small, and it is likely that G,(/, k) cannot be determined as 

accurately as in conditions of high SNR. In this situation, it is not so advantageous 
10 to amplify the Wiener filter output and the estimated Wiener filter should be 
maintained in the form it was originally estimated in step 3. 

To take into account these two contradictory requirements that exist in different 
SNR conditions, the Wiener filter determined in step 3 is modified according to: 

15 

G a (f,k) = G l (f,k) Min[Khini - G * {f > k)] (26) 

to produce a Wiener filter G a (f,k) to be used in estimation of the final output. 
G a (f,k) is a function of G x (f,k) . 

20 

Equation 26 exploits the fact that a function such as y = x l ~ x (x>0) provides 
amplification when x is less than one. It therefore fulfils the requirement of 
providing more amplification in good SNR conditions and less amplification in 
conditions of low SNR. 

25 

The variable Kb(f) can take values between 0 and 1 and is included in the 
exponent of Equation 26 in order to enable the use of different (e.g. 
predetermined) amplification levels for different frequency bands f, if desired. 

30 



Step 7: 



Selection of the Level of Noise Reduction 



Si 



pent 
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In this step, the desired level of noise reduction is selected. For the Wiener filter 
given in Equation 11, the corresponding ideal temporal output has the form 
j(0 = s(0 + Recalling that the noisy input signal has the form 

5 x(t) = s(t) + n(t) , the noise reduction provided by the filter is theoretically about 
20*log[£] dB. This result can be justified by considering the ratio of the noise level 
in the input signal to that in the output signal (i.e. the signal obtained after noise 
suppression). This ratio is simply £-n(0 / n(t) , which, when expressed as a 
power ratio in decibels, becomes 20 1og[£] dB. Consequently, the factor 0<£<1 
10 corresponds to the noise reduction introduced by the filter. 

Having chosen a desired noise reduction level and determined the value of £ 
necessary to achieve that noise reduction (e.g. for -12 dB noise reduction, 
£ = 0.25), a factor ?j is determined such that: 



G, (/, *) + 7/ • (1 - G 1 (/,*)) <=> p s(f> k ) + ^ p n(f> k "> . (27) 

/>,(/,*)+ P B (/,*) 



Equation 27 presents a way of relating a Wiener filter optimised to provide an 
output that includes only masked noise to a Wiener filter that provides an output 

20 including a certain amount of permitted noise. According to steps 1 - 3, the Wiener 
filter G, (f,k) is constructed so as to provide an estimate of the speech component 
of a noisy speech signal plus an amount of noise which is effectively masked by 
the speech component. Thus, in the condition where a certain amount of noise is 
permitted (desired) in the output, the Wiener filter must be modified accordingly. In 

25 Equation 27, G { (f,k) represents the Wiener filter optimised in step 3 to provide an 

P (f k) + £ • P (f k) 

output that contains speech-masked noise. The term * w — - — i^Lj^J. 

P s (f,k) + P n (f,k) 

represents a Wiener filter that provides an amount of noise reduction which 
produces an output signal containing speech and a desired/permitted amount of 
noise. The term 77 -(l-G, (/,&)) thus represents an amount of non-masked noise 
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and is essentially the difference between v W ' — - — ±i±J-J. and G. (/,£). Taking 

P s (f,k) + P n (f,k) 

into account the fact that G x {f,k) contains noise at a level of about (1-or) times 
the noise present in the original noisy speech signal, the following relationship 
between a , 77 and £ is tr ue: 



l-a + r/a&Z (28) 



Step 8: Estimation of the Final Estimated Wiener Filter 



10 Using Equations 16, 26 and 28, the final Wiener filter G(f,k) to be applied to the 
input is given by: 



a 



ifa>{\-£) 77- 
else 7] = 0 

[G(f,k) = G a (f,k) + 77 • (1 - G, (/,*)) 



(29) 



1 5 Although 77 depends on a, and has a different value for each frequency bin / of 
each frame k , the overall noise reduction level is maintained constant around 
20-log[£] dB. 



Alternatively, steps 1 to 8 could be implemented using formulae involving signal- 
20 to-noise ratio formulas. In the detailed implementation of steps 1-8, presented 
above, the discussion was based on calculations of noise psd functions, speech 
periodograms and input power (periodogram + psd). However, an alternative 
representation can be obtained by dividing Equation 11 and/or Equation 13 by the 
noise psd. This alternative representation requires estimation of a (signal+masked 
25 noise)-to-noise ratio, instead of a speech periodogram. 



An algorithm 50 embodying the invention is shown in Figure 5. The algorithm 50 is 
shown divided into a set of steps 52 which are an adaptive process and a set of 



+ 
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steps 54 which are a non-adaptive process. The adaptive process uses a 
computation of the Wiener filter to re-compute the Wiener filter. Accordingly, the 
step of the computation of the Wiener filter is common both to the adaptive 
process and to the non-adaptive process. 

5 

This Wiener filter calculation is also suitable for minimising the residual echo in a 
combined acoustic echo and noise control system including one sensor and one 
loudspeaker. 

10 While preferred embodiments of the invention have been shown and described, it 
will be understood that such embodiments are described by way of example only. 
For example, although the invention is described in a noise suppressor located in 
the up-link path of a mobile terminal, that is providing noise suppressed signal to a 
speech encoder, it can equally be present in a noise suppressor in the down-link 

15 path of a mobile terminal instead of or in addition to the noise suppressor in the 
up-link path. In this case it could be acting on a signal being provided by a speech 
decoder. Furthermore, although the invention is described in a mobile terminal, it 
can alternatively be present in a noise suppressor in a communications network 
whether used in relation to a speech encoder or a speech decoder. 

20 

Numerous variations, changes and substitutions will occur to those skilled in the 
art without departing from the scope of the present invention. Accordingly, it is 
intended that the following claims cover all such equivalents or variations as fall 
within the spirit and scope of the invention. 



